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[l] As an airless body in space with no global magnetic 
field, the Moon is exposed to both solar ultraviolet radiation 
and ambient plasmas. Photoemission from solar UV radia- 
tion and collection of ambient plasma are typically opposing 
charging currents and simple charging current balance pre- 
dicts that the lunar dayside surface should charge positively; 
however, the two ARTEMIS probes have observed energy- 
dependent loss cones and high-energy, surface-originating 
electron beams above the dayside lunar surface for extended 
periods in the magnetosphere, which are indicative of 
negative surface potentials. In this paper, we compare 
observations by the ARTEMIS PI spacecraft with a one- 
dimensional particle-in-cell simulation and show that the 
energy-dependent loss cones and electron beams are due 
to the presence of stable, non-monotonic, negative poten- 
tials above the lunar surface. The simulations also show that 
while the magnitude of the non-monotonic potential is 
mainly driven by the incoming electron temperature, the 
incoming ion temperature can alter this magnitude, espe- 
cially for periods in the plasma sheet when the ion temper- 
ature is more than twenty times the electron temperature. 
Finally, we note several other plasma phenomena associated 
with these non-monotonic potentials, such as broadband 
electrostatic noise and electron cyclotron harmonic emis- 
sions, and offer possible generation mechanisms for these 
phenomena. Citation: Poppe, A. R., J. S. Halekas, G. T. Delory, 
W. M. Farrell, V. Angelopoulos, J. P. McFadden, J. W. Bonnell, 
and R. E. Ergun (2012), A comparison of ARTEMIS observations 
and particle-in-cell modeling of the lunar photoelectron sheath in 
the terrestrial magnetotail, Geophys. Res. Lett., 39, L01102, 
doi: 10.1 029/20 11GL050321. 

1. Introduction 

[ 2 ] The surface of the Moon is exposed to a variety of 
ambient plasma, and in combination with solar ultraviolet- 
induced photoemission and the presence of surface remanent 
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magnetic fields, determines the plasma environment imme- 
diately above the surface. Early theoretical work on the 
charging of exposed surfaces in space, including the Moon, 
Mercury and asteroids, quickly recognized that for some 
combinations of photoelectron and ambient plasma para- 
meters, multiple solutions existed for the electrostatic 
potential structure above the surface [ Guernsey and Fu, 
1970; Fu, 1971; Grard, 1997; Nitter et al., 1998]. These 
solutions were classified as being either monotonic, where 
the potential either steadily increases or decreases from 
infinity to the surface, or non-monotonic, where the potential 
is increasing in one region and decreasing in another. 
Simulations of photoelectric charging of spacecraft have 
indicated the likely presence of non-monotonic potentials in 
situations with large photoemission and a relatively tenuous 
ambient plasma [Ergun et al., 2010]. These works also 
suggested, somewhat un-intuitively, that the non-monotonic 
potential may in fact be energetically preferred over the 
monotonic solution. In-situ measurements by the Lunar 
Prospector (LP) Electron Reflectometer (ER), which orbited 
the Moon and could sense the lunar surface potential via 
electron reflectometry, suggested that non-monotonic poten- 
tials were indeed present [Halekas et al., 2008b, 2012]. By 
comparing with a one-dimensional particle-in-cell simulation, 
these measurements were confirmed as resulting from non- 
monotonic potentials [Poppe and Horanyi, 2010; Poppe et al., 
2011 ]. 

[ 3 ] In this paper, we present a comparison of ARTEMIS 
PI measurements of the lunar dayside plasma environment 
with one-dimensional particle-in-cell simulations of the 
lunar photoelectron sheath for an extended period in the ter- 
restrial magnetotail. By making use of ARTEMIS’ extensive 
suite of particles and fields instruments [Angelopoulos, 
2012], including ion measurements (which were absent on 
LP), we are able to improve our particle-in-cell model and 
show excellent agreement with the ARTEMIS observations. 
This agreement significantly improves our understanding of 
the nature and variability of non-monotonic potentials and 
points to their likely ubiquity at airless bodies throughout the 
solar system. Additionally, we show that the presence of 
photoelectron beams associated with non-monotonic poten- 
tials generates large amounts of broadband electrostatic 
noise, most likely due to streaming instabilities. 

2. ARTEMIS Observations 

[ 4 ] The ARTEMIS mission consists of two identical 
probes in elliptical orbits around the Moon, with PI orbiting 
prograde and P2 orbiting retrograde [Angelopoulos, 2012]. 
In this study, we focus on measurements by PI on July 16, 
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2011, of low-energy electrons and ions, and magnetic and 
electric field wave spectra, using the Electrostatic Analyzer 
(ESA) [ McFadden et al ., 2008], Electric Field Instrument 
(EFI) [Bonnell et al , 2008] and Search Coil Magnetometer 
(SCM) [Roux et al , 2008]. At this time, the Moon was 
located deep in the terrestrial magnetotail at GSE coordi- 
nates of approximately, [—59.5, —13.5, 3.80] R E and the PI 
orbit placed the spacecraft above the dayside lunar surface 
throughout the entire measurement period. 

[ 5 ] Figure 1 shows an overview of PI observations 
between 11:50 and 12:25 UTC, when the spacecraft was 
magnetically connected to the dayside lunar surface, assum- 
ing a straight-line extrapolation. While the incident energy 
flux is typical of the terrestrial magnetotail, the flux coming 
from the lunar surface is drastically different, showing both a 
relative deficit of high-energy electrons, which are lost to the 
lunar surface, and a cold beam of electrons accelerated up 
from the lunar surface. The beam energy varies throughout 
this time, including some periods in which the beam dis- 
appears, such as 1 1:59 to 12:01. In the magnetic and electric 
field FFT wave spectra, a host of electrostatic and electro- 
magnetic waves are seen. From 11:50 to 11:58, and again 
from 12:01 to 12:14, both spectra show the presence of a 
narrow-band electromagnetic wave at frequencies varying 
between 30 and 200 Hz. These waves are most likely elec- 
tromagnetic whistler waves, unrelated to the electrostatic 
potentials, and thus, we defer discussion of these to future 
work. The electric field FFT wave spectra show the presence 
of broadband electrostatic noise (BEN) throughout this time 
period, correlated with the presence of the up-going electron 
beam. Additionally, between approximately 1 1:57 and 12:00 
there are a series of narrowband electrostatic emissions at 
approximate harmonics of the electron cyclotron frequency. 

[6] Figures 2a and 2b show electron pitch angle-energy 
spectrograms in units of differential energy flux for two 
periods of distinctly different observations. Figure 2a shows 
the spectrogram at 11:59:39, with down-going plasma tem- 
peratures of T e = 59 eV and T L = 762 eV for electrons and 
ions, respectively, while Figure 2b shows the spectrogram at 
12:04:04, with T e = 247 eV and T t = 1,310 eV. In both 
cases, the down-going plasma sheet electrons are isotropic, 
even at differing temperatures; however, the up-going 
spectra show distinct anisotropies as a result of combined 
electrostatic and magnetic reflection near the lunar surface 
[Halekas et al , 2002, 2008b]. In both Figures 2a and 2b, an 
energy-dependent loss cone is seen for high-energy, field- 
aligned electrons that are absorbed by the lunar surface after 
penetrating both the lunar crustal remanent magnetic fields 
and the electrostatic potential. Additionally, both observa- 
tions show up-going fluxes greater than the corresponding 
down-going flux at certain pitch angles and energies above 
the spacecraft potential (denoted as a vertical white line in 


Figures 2a and 2b). Below this potential, spacecraft- 
generated photoelectrons contaminate the detector [Halekas 
et al , 2005]. In Figure 2a, the surface-originating flux 
manifests as an increase in up-going fluxes in the lowest 
non-contaminated bin centered at approximately 60 eV 
across 0-90° pitch angle, while in Figure 2b, the surface- 
originating flux is a cold beam at energies between 150— 
200 eV and pitch angles less than 30°. Spectra similar to 
Figure 2b have been previously shown to be indicative of 
negative, non-monotonic lunar surface potentials with 
respect to the ambient plasma [Halekas et al , 2005, 2008b, 
2011; Poppe et al, 201 1]. Figure 2a, in which the up-going 
photoelectrons do not develop into a focused beam in pitch 
angle, is still indicative of non-monotonic potentials, but at 
a much smaller magnitude than Figure 2b, assuming that the 
lunar photo electrons are emitted as a roughly-Maxwellian 
distribution with a temperature of approximately 2.2 eV 
[Feuerbacher et al, 1972]. 

3. Model / Data Comparisons 

3.1. PIC Model Description and Results 

[ 7 ] Previous comparisons between LP ER data and a one- 
dimensional particle-in-cell model of dayside negative sur- 
face potentials in the terrestrial plasma sheet have concluded 
that observations similar to Figures 1 and 2 are most likely 
due to the presence of stable, non-monotonic potential layers 
above the lunar surface [Poppe et al , 2011; Halekas et al , 
2011]. In that study, only the electron parameters were 
explicitly measured, as LP lacked any ion instrumentation. 
With the use of ARTEMIS data, our model, as described by 
Poppe and Horanyi [2010] and Poppe et al [201 1], can now 
be improved to include both the electron and ion tempera- 
tures, which are significantly different during this time 
period. Shown in Figure 3 is the PIC-modeled electrostatic 
potential as a function of height above the lunar surface for 
three values of the ion temperature, 7) = 500, 1000 and 
2000 eV, with a constant electron temperature of T e = 500 eV 
and density at infinity of n ps = 10 5 m -3 . For all three cases, 
the potential is non-monotonic, yet with increasing ion tem- 
perature, the magnitude of the potential drop from infinity 
decreases as the ions provide more current to the surface. 
From the model results, we can derive an expression for the 
photoelectron beam energy, E b , equivalent to the absolute 
value of the electrostatic potential drop from infinity to the 
minimum, as a function of the incoming ion and electron 
temperatures, T e and T h given by, 

E b (T e , Ti) = 0 . 7127 ; - 0 . 0537 ; + 5 . 92 , ( 1 ) 

where all quantities are in units of eV. The incoming plasma 
sheet density, n ps , does play a role in modulating the 


Figure 1. ARTEMIS observations from a magnetically connected period on July 1 6, 20 1 1 , in the terrestrial magnetosphere. 
Shown are the differential energy flux for electrons coming up from (0—30° pitch angle) and going down to (90—180° pitch 
angle) the lunar surface, respectively, measured in eV/cm 2 /s/str/eV, the downgoing ion and electron temperatures, the high- 
frequency FFT wave spectra for the magnetic and electric fields, respectively, the low frequency magnetic wavelet, the dis- 
tance to the lunar surface along the magnetic field line in lunar radii, and the solar zenith angle (SZA), azimuthal (phi) and 
polar (theta) angles of the connected magnetic field line, assuming a straight-line extrapolation. Additionally, over-plotted in 
white on the differential electron energy fluxes (top two panels) is the spacecraft potential, over-plotted in red on the FFT 
spectra is the electron cyclotron frequency, and over-plotted in red on the wavelet is the proton cyclotron frequency. The 
two dashed lines denote individual observations shown in Figure 2. 
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electrostatic potential; however, the effect is negligible for 
these observations as the density varied within a factor of 
two, compared to the temperatures, which have variations of 
approximately an order of magnitude. 


3.2. Comparison to ARTEMIS Data 

[8] To explicitly compare the model predictions with the 
ARTEMIS observations, we use equation (1) with the 
ARTEMIS-measured T e and T t as a function of time to yield 
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Figure 2. Electron pitch angle-energy spectrograms shown 
in units of differential electron energy flux for two time 
periods, (a) 1 1 :59:39 and (b) 12:04:04. The spacecraft poten- 
tial for each case is plotted as the vertical white line. 

a PIC model-predicted beam energy as a function of time. 
Figure 4 (top) shows the ratio of up-going to down-going 
electron differential energy flux with the PIC model- 
predicted beam energy (summed with the spacecraft poten- 
tial), for two conditions: (1) T t = T e (red line), and (2) T t 
equivalent to the true ion temperature, as measured by 
ARTEMIS (black line). The PIC model-predicted beam 
energies for both conditions show excellent agreement with 
the observed beam, especially for times between 12:02 and 
12:25, where the ion temperature is only two to five times 
the electron temperature. For times earlier than 11:59, the 
observed beam appears at energies slightly less than the 
predicted value. During this time, there are indications of 
electromagnetic wave activity near the proton cyclotron 
frequency, which is indicative of ion dynamics, perhaps due 
to ions of lunar origin; however, the Moon was in the mag- 
netotail and no significant convection electric field exists to 
pick up any ions generated near the lunar surface (via photo- 
ionization or sputtering, for example) and accelerate them to 
measurable energies. Nonetheless, one plausible mechanism 
that could couple the presence of multiple ion species and the 
cold photoelectron beam is the generation of electromagnetic 
ion cyclotron waves, as described by Temerin and Lysak 
[1984]. Indeed, the observed deviations of the photoelec- 
tron beam energy from the model may be telling us about 
near-surface particle populations and dynamics. 

[9] While the photoelectron beam is a prominent effect 
throughout most of this time period, it does disappear, in 
agreement with the PIC model, from approximately 11:5 9— 
12:02. Concurrent with this time period, the plasma sheet 
electron temperature (third panel of Figure 1) decreases 
smoothly from approximately 200 eV to a minimum around 
50 eV before smoothly increasing back to its previous value. 


The plasma sheet ion temperature remains relatively con- 
stant during this time at approximately 1000 eV. By setting 
the left-hand side of equation (1) equal to zero, we can 
determine values of T e and T t for which the predicted beam 
energy is less than or equal to zero (implying that the beam 
vanishes). By first defining the temperature ratio, R = TJTe , , 
we solve equation (1) for a somewhat implicit, model- 
predicted condition for the beam to vanish, given by, 

R > 13.4+ III.7T (2) 

1 e 

[10] Figure 4 (middle) shows the temperature ratio as a 
function of time, with a critical ratio, defined by 
R c = 13.4 + Hl.777 1 , over plotted in red. During most of 
this time, the temperature ratio is below the critical value and 
the beam is seen in the up-going energy fluxes. Commen- 
surate with the decrease of the incoming electron tempera- 
ture, the ratio increases above R c and correspondingly, both 
the model and the data show the beam disappear, although, 
there does exist a small amount of un-focused, upward 
energy flux in the lowest energy bins above the spacecraft 
potential that serves to balance the plasma sheet current to the 
lunar surface. During this time, the plasma sheet electron 
current to the lunar surface decreases as the electron temper- 
ature decreases and thus, the net current to the lunar surface 
increases. The non-monotonic potential correspondingly 
decreases in magnitude to reduce the incoming ion current and 
maintain current balance at the surface, to the extent that the 
photoelectrons escaping from the lunar surface no longer 
experience significant acceleration and pitch-angle focusing as 
they travel up to the spacecraft (see Figure 2a). 

[11] Finally, Figure 4 (bottom) shows the electric field 
wave power in the 732 Hz bin, which, based on Figure 1 , can 
be used as a proxy for the presence of BEN. The BEN is 
present at times when the beam is at relatively high energies, 
yet ceases when the beam has either disappeared (11:5 9— 
1 2 : 0 1 ) or when the beam is at a very low energy level (12:13- 
12:25). One possible mechanism for the generation of both 
the BEN and the previously mentioned electron cyclotron 
harmonic emissions may be the coupling of the electron 
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Figure 3. The potential as a function of height above the 
lunar surface from the PIC model for three different cases 
of plasma sheet ion temperature, T t = 500, 1000, 2000 eV, 
for constant electron temperature, T e = 500 eV, and density, 
n ps = 10 5 m -3 . 
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2011 Jul 16 

Figure 4. (top) The up-going to down-going differential electron energy flux ratio with the model-predicted beam energy 
with T t = T e and the true T t taken from ARTEMIS measurements (red and black lines, respectively), and the spacecraft poten- 
tial (white line) over-plotted, (middle) the temperature ratio ( TJT e ) with the critical ratio, R c , as defined in the text (red line), 
and (bottom) the electric field wave power in the 732 Hz bandpass for the same time period as Figure 1. 


acoustic instability from the streaming photo electrons and 
the instability of the loss cone anisotropy, similar to that 
outlined by Roeder et al. [1991], especially given the high 
degree of correlation between the BEN and the presence of 
the photoelectron beam. We suggest that the presence of both 
BEN and electron cyclotron harmonic emissions near airless 
bodies can be used as a diagnostic in searching for surface- 
generated electron beams and possible non-monotonic 
potentials at other airless bodies [ Gurnett et al . , 1981; Kurth 
et al., 1987, 2001; Santolik et al., 2011]. 

4. Summary 

[ 12 ] We have presented measurements by ARTEMIS PI 
of energy-dependent loss cones and cold, up-going field- 
aligned beams of photoelectrons, both of which indicate that 
the dayside lunar surface is at a negative potential with 
respect to the spacecraft. Using a one-dimensional particle- 
in-cell code, we have shown that these measurements are 
due to stable, non-monotonic potentials that arise due to the 
interaction of the lunar photoelectron and terrestrial plasma 
sheaths. The extended time period and plasma temperature 
range over which the PIC model agrees with the ARTEMIS 
measurements significantly improves our understanding of 
non-monotonic potentials above airless bodies, which has 
previously only been compared for a single Lunar Prospector 


measurement [Poppe et al., 2011]. Indeed, the wide tem- 
perature range over which the PIC model has now been 
validated against in-situ observations, from 10 eV in the 
solar wind [Poppe and Horanyi, 2010; Halekas et al., 201 1] 
to over 1 keV in the terrestrial plasma sheet [Poppe et al., 
2011], implies that non-monotonic potentials may exist at 
nearly all airless, non-magnetized bodies throughout the 
solar system. The presence of a multitude of electrostatic and 
electromagnetic waves associated with the non-monotonic 
potentials, loss cones, and photoelectron beam further 
enriches this environment, and provides additional means of 
searching for non-monotonic potentials throughout the solar 
system. We suggest that these phenomena are not unique to 
the Moon and should occur at many airless bodies, includ- 
ing Mercury, asteroids, and the moons of the giant planets 
[Grard, 1997; Roussos et al., 2010; Santolik et al., 2011]. 
The non-monotonic potentials also have implications for 
other lunar-related phenomena, including the electrostatic 
levitation of micron-sized dust grains [Poppe and Horanyi, 
2010; Griin et al., 2011] and the interaction of plasma with 
lunar crustal remanent magnetic fields [Harnett and Winglee , 
2003; Halekas et al, 2008a]. Finally, we also suggest that 
the continuous presence of such non-monotonic potentials, 
and the complex electric fields that these potentials give rise 
to, may have significant effects on the trajectories and 
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densities of sputtered and photo-ionized pick-up ions above 
the day side lunar surface. 
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